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(57) Abstract 

In accordance with the method a turbid medium (8) is scanned by displacing an opdcal signal source and a coiresponding optical 
detector (5) from one respective spatial location to another about the turbid medium, each spatial location being associated with conesponding 
input regions on opposite faces of the turbid medium (8). The optical signal source directs an optical signal (1) to each of said input regions. 
The optical detector (5) sends corresponding detector output signal data to an image processing means (5a) for generating a set of nnage 
data of the turbid medium (8). The improvement comprises generating in addition to a primary set of image data, secondary and/or tertiary 
image data which are derived by scanning the tuibid medium (8) using input regions and output regions of different sizes. The various sets 
of image data obtained may be manipulated to highlight the differences or similarities between the images (e.g. to obtain a composite set 
of image data indicative of an anomaly(ies) in the body of the turbid medium (8). 
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DEPTH DISCRIMINATION OF HETOROGENEITIES IN TURBID MEDIA 



The present invention relates to optical imaging of turbid media, e.g. breast tissue. 

Optical maounography is likely to become a relevant technology for breast cancer detection. 
5 Light scattering is the main problem that limits the potential of efficiently detecting small 

lesions especially when located in the central region of the breast Structures close to surface 
are easily detected by optical techniques while more deep structures are more difficultly 
detected. A technique to preferentially probe in the central region of the breast is thus of 
interest. 

10 The present text will refer by way of example only to breast tissue; however, other types of 
tuibid media (i.e. other types of (living) tissue) may also be imaged using the techniques 
described herein. 

Optical imaging of turbid media has been the subject of much research activity and has seen 
an increase in interest since the early 1990*s. This type of imaging is based on the fact that the 

15 propagation of light in a turbid medium depends on the absorption and scattering properties 
of the medium. Absorption results from energy level transitions of the constituent atoms and 
molecules in the medium. The absorption property of the medium is quantified by its 
absorption coefficient m defined as the probability of a photon being absorbed per 
infinitesimal pathlength. Scattering results from variations in the index of refi-action of the 

20 different structures present in the mediiun. In a highly diffusive medium, scattering is 

quantified by the reduced scattering coefficient defined as the probability of a photon 
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being isotropically scattered per infinitesimal pathlength. Characteristics such as intensity, 
coherence and polarization of the incident light change as it is absorbed and scattered by the 
medium resulting in diffuse transmittance of the light. In particular, scattering causes a 
collimated laser beam to spread over a sizeable volume element. This complicates the 

S imaging of a turbid medium. Special imaging modalities must be implemented to offset the 
detrimental light dif&sion. For example, time-resolved methods use short light pulses from a 
fast laser source to illiuninate the mediiun. The emergent light is collected by a fast detector 
capable of reproducing its time variation, which can provide further information about the 
turbid medium. A simple data processing approach in this case is time-gating, by which only 

10 the earliest part of the output light pulses is used to produce an image. This amoimts to using 
only the light with the straightest trajectory through the scattering medium, thus improving 
spatial resolution (please see J. C. Hebden and R. A. Kruger, "Transillumination imaging 
performance: Spatial resolution simulation studies", Med. Phys. 17, 41^7 (1990)) Reference 
may also be made to the following documents for additional information with respect to 

IS optical scanning techniques: S. B. Colak, D. G. Papaioannou, G. W. Hooft, M. B. Van der 
Mark, H. Schomberg, J. C. J. Paasschens, J. B. M. Melissen, and N. A. A. J. Van Asten, 
'Tomographic image reconstmction from optical projections in light-diffusing media", Appl. 
Opt. 36, 180-213 (1997); M. S. Patterson, B. Chance, and B. C. Wilson, "Time-resolved 
reflectance and transmittance for the noninvasive measurement of tissue optical properties", 

20 Appl. Opt. 28, 2331-2336 (1989); D. Contini, F. Martelli. and G. Zaccanti, "Photon 

migration through a turbid slab described by a model based on diffusion approximation. I. 
Theory", Appl. Opt. 36, 4587-4599 (1997); M. Morin, S. Chatigny, A. Mailloux, Y. 
Painchaud, and P. Beaudry, "Time-domain perturbation analysis of a scattering slab", in 
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Optical Tomography and Spectroscopy of Tissue IIIj B. Chance, R. R. Alfano, and B. J. 
Tromberg, eds., Proc. SPIE 3597, 67-78 (1999); Y. Painchaud, A. Mailloux, M. Morin. S. 
Verreaull and P. Beaudry, "Time-domain optical imaging: discrinlination between scattering 
and absorption", AppL Opt. 38, 3686-3693 (1999); as well as U.S patent number 5,808,304). 

S The strong interest in optical imaging of scattering media stems from the need for biomedical 
diagnostic techniques that are safe and non-invasive. The optical properties of biological 
tissues are at the heart of optically based biomedical diagnostic techniques. As for the general 
case of a turbid medium, the manner in which light propagates through tissue depends on its 
absorption and scattering properties. Thus, if abnormal tissue can be said to differ from 
10 normal in its absorption or scattering of light for some physiological or morphological reason,, 
it then becomes possible to optically differentiate between normal and abnormal conditions. 
, A specific application is optical mammography where tumors could be differentiated from 
normal breast tissue on the basis of optical properties. 

It would, for example, be advantageous to have a mechanism for the detection of an 
1 5 anomaly(ies) in turbid medium such as for example a tumour in tissue such as breast tissue. 

Accordingly, the present invention in one aspect provides in a method for scanning a turbid 
medium for generating an image thereof for the detection of one or more anomalies contained 
within the turbid medium, said turbid medium having a first face and an opposite second face, 
wherein said turbid medium is scanned by displacing an optical signal source and a 
20 corresponding optical detector from one respective spatial location to another about the 

turbid medium, each spatial location being associated with a corresponding input region on 
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said first face and a corresponding output region on the opposite second face, said optical 
signal source directing an optical signal (e.g. laser beam) to each of said input regions, said 
optical detector in response to optical signals detected from each of said output regions 
sending corresponding detector output signal data to an image processing means for 
5 generating a set of image data of the turbid medium, the improvement comprising 

generating a primary set of image data derived by scanning said turbid medium using 
for each input region a first input region of predetermined size and using for each output 
region a first output region of predetermined size, and 

generating a secondary set of image data derived by scanning said turbid medium 
10 using for each input region said first input region of predetermined size and using for each 
output 

The present invention in another aspect provides in a method for scanning a turbid medium 
for generating an image thereof for the detection of one or more anomalies contained within 
the turbid mediiun, said turbid medium having a first face and an opposite second face, 

15 wherein said turbid medium is scanned by displacing an optical signal source and a 

corresponding optical detector from one respective spatial location to another about the 
turbid medium, each spatial location being associated with a corresponding input region on 
said first face and a corresponding output region on the opposite second face, said optical 
signal source directing an optical signal (e.g. laser beam) to each of said input regions, said 

20 optical detector in response to optical signals detected from each of said output regions 
sending corresponding detector output signal data to an image processing means for 
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generating a set of image data of the turbid medium, the improvement comprising 

generating a primary set of image data derived by scanning said turbid medium using 
for each input region a first input region of predetermined size arid using for each output 
region a first output region of predetermined size, and 

5 generating a further set of image data derived by scanning said turbid mediiun using 

for each input region a second input region of predetermined size greater than that of said first 
input region and using for each output region said first ou^ut region of predetermined size. 

The present invention in another aspect provides in a method for scanning a turbid medium 
for generating an image thereof for the detection of one or more anomalies contained within 

10 the turbid medium, said turbid medium having a first face and an opposite second face, 
wherein said turbid medium is scanned by displacing an optical signal source and a 
corresponding optical detector fi-om one respective spatial location to another about the 
turbid medium, each spatial location being associated with a corresponding input region on 
said furst face and a con-esponding output region on the opposite second face, said optical 

15 signal source directing an optical signal (e.g. laser beam) to each of said input regions, said 
optical detector in response to optical signals detected fi^om each of said output regions 
sending corresponding detector output signal data to an image processing means for 
generating a set of image data of the turbid medium, the improvement comprising 

generating a primary set of image data 1^ derived by scanning said turbid medium 
20 using for each input region a first input region of predetermined size and using for each 
output region a first output region of predetermined size, 

'5 
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generating a secondary set of image data Ib derived by scanning said turbid medium 
using for each input region said first input region of predetermined size and using for each 
output region a second output region of predetermined size greater than that of said first 
output region and 

5 generating a tertiary set of image data Ic derived by scanning said tuibid medium 

using for each input region a second input region of predetermined size greater than that of 
said first input region and using for each output region said first output region of 
predetermined size. 

In accordance with the present invention the various sets of image data may be generated in 
10 any desired or necessary order. 

In accordance with the present invention, the various sets of image data obtained as described 
herein may be manipulated in any manner (i.e. subjected to a data processing technique) so as 
to for example highlight the differences and/or similarities between the images (e.g.. to obtain 
a composite set of image data indicative of an anomaly(ies) (e.g. breast tumour) in the body 

IS of the turbid medium); such manipulation or data processing treatment may be carried out 

by simple subtraction or by exploiting more sophisticated techniques such as for example K- 
L transforms. Thus, for example, the method(s) described herein may comprise (a) 
combining said primary set of image data and said secondary set of image data so as to 
generate a product set of image data wherein the similarities and/or differences between the 

20 said primary set of image data and said secondary set of image data are highlighted.; (b) 

combining said primary set of image data and said further set of image data so as to generate 

6 
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a product set of image data wherein the similarities and/or differences between the said 
primary set of image data and said further set of image data are highlighted.; (c) combining 
said primary set of image data said secondary set of image data and said tertiary set of image 
data so as to generate a product set of image data wherein the similarities and/pr differences 
5 between said primary set of image data, said secondary set of image data and said tertiary set 
of image data are highlighted; etc.. Thus, more particularly, for example, a composite set of 
image data may be generated by subtracting the secondary set of image data from the primary 
set of image data; a composite set of image data may be generated by subtracting the further 
set of image data from the primary set of image data; a composite set of image data may be 
1 0 generated in accordance with the equation Ib Ic - 

In accordance with the present invention a system for performing optical imaging 
measurements on a turbid mediimi in accordance with the method(s) as described herein may, 
for example, comprise 

any suitable (known) means for generating a coUimated optical signal (e.g. laser) 

1 5 directing means for directing said optical signal to an input area or region of 

predetermined size on a face on one side of the turbid medium (e.g. using a suitable optical 
fiber) 

detection means for detecting a optical output signal from a second area or region of 
predetermined size on the face of the turbid mediiun on the opposite side of said turbid 
20 medium (e.g. photomuliplier tube, streak camera, etc.); the detecrion means may include 

optical fibers for directing the optical signal to the component of the detection means which 

7 
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may convert the optical signal to an output (e.g. analogue or digital electrical) signal; and 

means (e.g. computer) for processing said optical output signals (in any known 
manner) so as to obtain a set of image data (e.g. data which may be projected onto the screen 
of a computer monitor, be printed on a printer, be further processed, stored in computer 
5 memory, etc.). 

The system may also of course incorporate any suitable mechanism or means for 
manipulating the size of the input and ou^ut regions as referred to herein. 

It is to be understood herein (i.e. for the purposes of the present invention) that the optical 
signal(s) may be of electromagnetic radiation which is not only in the visible but as well as in 
10 the infra-red or the near infra-red part of the spectmm keeping in mind the purpose of the 
signal; the initial optical signal may for example be a laser beam. 

The technique of the present invention may as mentioned above involve two or three 
(successive) scans of the turbid medium (e.g. breast tissue) so as to provide two or three 
projection images. These two or three images may be used to image stmctures according to 
15 their depth inside the breast. The first image shows all the stmctures as in the usual optical 
techniques. A second image preferentially highlights the stmctures near the input surface. 
Another, third image preferentially highlights the structures near the output surface. By a 
combination of two or three of the images, it is possible to highlight the stmctures in the 
central region of the breast and attenuate the stmctures near the surfaces. 

20 In the figures which illustrate example embodiments of the present invention: 

8 
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Figure 1 illustrates in schematic fashion a Scanning system for imaging through a turbid 
medium. 

Figure 2 illustrates in schematic fashion a Scanning system for imaging through a turbid 
medium having a point-source/point-detector configuration 

Figure 2A shows a graph for the configuration of Figure 2 wherein the qualitatively indicates 
the transmitted light versus transverse position for each configuration 

Figure 3 illustrates in schematic fashion a Scanning system for imaging through a turbid 
medium having a point-source/large-area detector configuration 

Figure 3 A shows a graph for the configuration of Figure 3 wherein the qualitatively indicates 
the transmitted hght versus transverse position for each configuration 

Figure 4 illustrates in schematic fashion a Scanning system for imaging through a turbid 
medium having a large-area source/point-detector configuration 

Figure 4A shows a graph for the configuration of Figure 4 wherein the qualitatively indicates 
the transmitted light versus transverse position for each configuration; 

Figure 5 illustrates a firont view of a phantom; 

Figure 6 illustrates a side view of the phantom of figure 5; 

Figure 7 is a graphic illusu-ation of Spatial profile combinations used to enhance longitudinal 

9 
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discrimination; 

Figure 8 schematically illustrates an optical fiber arrangement for a point-source/point 
detector configuration 

Figure 9 schematically illustrates an optical fiber arrangement for a point-source/large-area 
5 detector configuration 

Figure 10 schematically illustrates an optical fiber arrangement for a large-area source/point- 
detector configuration 

Figure 1 1 is an illustration of the spatial profile combination for Pa-Pb of a liquid phantom 

Figure 12 is an illustration of the spatial profile combination for Pa-Pc of a liquid phantom 

10 Figure 1 3 is an illustration of the spatial profile combination for Pb + Pc -Pa of a liquid 
phantom 

Figure 14 is an illustration of the spatial profile for Pa of a liquid phantom 

In biomedical optical imaging, two types of images can be generated: 3D reconstmcted 
images and 2D projection images. The foraier are produced using tomography which is 
15 typically based on a multi-point geometry involving a large number of detectors ( see S* B. 
Colak, D. G. Papaioannou, G. W. Hooft, M. B. Van der Mark, H. Schomberg, J. C. J. 
Paasschens, J. B. M. Melissen, and N. A. A. J. Van Asten, "Tomographic image 
reconstruction fi-om optical projections in light-diffusing media", Appl. Opt. 36, 180-213 

10 
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(1 997)). Its advantage is that 3D images are generated, however, measurements and 
reconstructions are potentially time-consuming. 2D projection images are based on the 
detection of the light that enters the medium through a small area, that propagates through the 
medium and that exits the medium over a small detection area facing the entrance area. The 
5 image is generated by measuring the transmitted light in such a manner at a large number of 
points as both the light beam incident on the medium and the detection device are scanned 
synchronously along the sur&ce of the medium. This imaging scheme is illusUiited in Figure 
1 wherein arrow 1 is indicative of a laser beam; wherein arrow 3 is indicative of the scanning 
direction of the beam 1 ; wherein the optical detector means is indicated by reference numeral 

10 5; wherein arrow 7 is indicative of the scanning direction of the optical detector means 5; and 
wherein the scattering medium is designated by the reference numeral 8. For the scanning 
system schematically illustrated in figure 1 (known) means are provided (not shown) whereby 
the laser beam 1 and the optical detector means 5 are respectively displaced in the direction of 
the arrows 3 and 7 in synchronise fashion, i.e. such that the optical detector 5 is always 

15 opposite the laser beam 1 . Although also not shown (in detail) in figure 1, the optical 

detector 5 (as well as the means for generating the laser beam and the means for directing the 
laser beam to the scan surface or face ) may be coupled in any (known) manner to a suitable 
(known) processing unit (e.g. computerised processing system) for ultimately generating or 
deriving an image (i.e. a set of image data) from the set of signal data obtained from the 

20 optical detector 5 during the scanning process. Figure 1 does show in block diagram 
fashion a processing unit 5a coupled to the optical detector 5 for processing the detector 
output signal data into image data. 

11 
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The scanning technique illustrated by the system in figure 1 has the advantage that it is fast 
and compatible with time-resolved measurements. However, information is limited to two 
dimensions, the detected light giving infomiation about a volume extending over the whole 
Une-of-sight joining the input point of the laser beam and the detector. This is illustrated in 
5 Figure 1, where the gray region 10 within the dotted lines 1 1 and 13 represents the volume 
through which detected photons have most likely propagated. The shape of this volume can 
be understood by considering that all photons enter the scattering medium at the same point 
and all detected photons leave it through a small area facing the detector. On the other hand, 
scattering allows detected photons to wander away from the direct line-of*sight joining the 

10 laser beam and the detector, this wandering being maximtun at the half-distance between the 
two. As illustrated in Figure 1, the scanning technique provides no longitudinal information, 
since the probed volume extends over the whole thickness of the scattering medium. 
Moreover, a small volume heterogeneity, with optical properties different from those of the 
scattering medium, is more likely to be visible when located close to the input point of the 

15 laser beam or the output point facing the detector. This is so because it is then probed by most 
detected photons, whereas a smaller fraction of detected photons propagate through the same 
heterogeneity standing in the middle of the scattering mediimi, where the probed volume is 
the widest. For optical mammography this can be troublesome, as a tumor located deep 
within breast tissue could be overshadowed by strongly absorbing stractures such as veins 

.20 located near the surface. 

The detection method in accordance with the present invention may provide for the basic 
evaluation of heterogeneity depth in a scattering medium (e.g. turbid medium) to be used in 

12 
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conjunction with optical imaging based on light transmittance measurement and the 
production of projection images. Depth assessment is necessary for a full 3-D localization of 
the heterogeneity. It can also aid visualization of stmctures located in the central plane of a 
turbid medium by identifying those structures located close to surfaces, which can then be 
taken out of images. Finally, knowledge of the longitudinal position of the heterogeneity can 
be used during data processing to improve the quantification of the optical contrast of the said 
heterogeneity with regards to the surrounding scattering medium. Within the context of 
manunography, this can help in tissue recognition. 

The present invention which may be referred to by way of example as Dual Spatial 
Integration (DSI) permits localization by separating the scattering medium into two or three 
longitudinal zones: a zone close to the input surface where the light is injected, a central zone 
and/or a zone close to the output surface where the light is detected^ By dividing the mediimi 
into two or three segments, stmctures close to surfaces can be extracted and set aside, and 
structures located in the bulk enhanced. 

Figure 2, 3 and 4 (as well as respective figures 2A, 3 A and 4A) illustrate (in schematic 
fashion) how detection surfaces and input beams of different sizes may be used to achieve a 
certain degree of longitudinal localization of an inclusion within a scattering slab, e.g. by 
appropriately modifying the input or output region sizes while maintaining the other 
variables (i.e. essentially the same) during the generation of each of the sets of image data 
during the scanning procedures. 



13 
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DSI may as mentioned above use three different source-detector configurations to achieve 
this sectioning; example configurations are illustrated schematically in figures 2, 3 and 4; for 
each of the shown configurations the same scauering medium is used, i.e. a scattering 
medium containing three inclusions at different depths (i.e. the three black dots or circles 
5 shown in these figures). 

Specifically, the first configuration (i.e. configuration 1) shown in figure 2 comprises a point 
source (i.e. small area input laser beam) 20 and a point detector (i.e. small area detector) 22 
which are used as the base configuration. The second configuration (i.e. configuration 2) 

10 shown in figure 3 comprises a point source (i.e. small area input laser beam) 20 and a large 

detector (i.e. large area detector) 24. The third configuration (i.e. configuration 3) shown in 
figure 4 comprises a large source (i.e. large area input laser beam) 26 and a point detector 
(i.e. small area detector) 22 are used as the base configiu'ation. For each of the configurations 
the laser beam is generated in any known suitable fashion; the light leaving the scattering 

15 medium is detected by any known suitable detector; the beam and detector (or alternatively, 
the medium itselO (is) displaced (i.e. so as to facilitate scanning of the medium) 
longitudinally in synchronise fashion in any (known) suitable fashion across opposed faces 
of the scattering medium; and images are generated from the detector signals by any suitable 
known processing system (e.g. computer system). For each of figures 2, 3 and 4: the gray- 

20 shaded rectangle is the turbid medium; the dark circles indicate the placement of the 

inclusions or heterogeneities in the turbid medium; the dark line above each rectangle refers 
to either a large area detector (long line) or a point detector (short line) - reference numbers 
22 and 24; the arrows indicate the type of source used: single arrow for a point source and 
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multiple arrows for a large area source, namely the elements specified by refereiice numbers 
20 and 26; and the thin lines within the turbid medium delimitate qualitatively the volume 
being probed by detected photons. The figures 2 A, 3 A and 4A show graphs wherein each 
respectively qualitatively indicates the transmitted light versus transverse position for each 
5 configuration. 

In accordance with the present invention the detection method may proceed as discussed 
hereinafter. The scattering medium is scanned with the scanning configuration shown in 
figure 2 and an image (1^) generated. In the image (Ia) so obtained, structures close to both 

10 the exit and entrance faces are highlighted relative to those within the bulk. . As the laser 
beam and detector are moved transviersely along the surface of the scattering mediuni, the 
inclusions close to the input and output surfaces have a stronger impact on the detected light, 
as discussed above. This is represented qualitatively by the graph shown in figure 2A, which 
depicts the detected light (signal) as a function of the position across the surface of the . 

15 scattering medium. A large area detector is then substituted for the point detector so as to 

obtain the second configuration shown in figure 3. This changes the volume that is probed by 
detected photons, as illustrated in Figure 3. The scattering medium is scanned with this 
second configuration and an image (Ib) generated. In this image (Iq), structures closest to the 
entrance face are highlighted, relative to those within the bulk or close to the exit face, 

20 because they are probed by a larger firaction of the detected photons. This is represented 
qualitatively by the graph shown in figure 3A, which depicts the detected light (signal) as a 
fimction of the position across the surface of the scattering medium. Thirdly, a large area 
source is substituted for the point source and a point detector replaces the large area detector 
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SO as to obtain the third configuration. This changes again the shape of the volume being 
probed by detected photons, as shown in Figure 4. The scattering medium is scanned with 
this third configuration and an image (Ic) generated. In this image (Ic), stractures closest to 
the exit face are now highlighted relative to those within the bulk or close to the entrance 
5 face. This is represented qualitatively by the graph shown in figure 4A, which depicts the 
detected light (signal) as a function of the position across the surface of the scattering 
medium. 

If images Ib and Ic (i.e. their image or signal data) are added and image 1^ (i.e. its image or 
10 signal data) subtracted, stmctures located within the bulk of the scattering medium can be 
more readily visualized than in images generated with the base configuration (Ia) alone* 

It is to be understood herein that the expression "small area" is refers to a surface whose 
diameter is relatively small compared to the scattering (turbid) medium thickness while the 
15 expression "large area" is to be understood as referring to a surface whose diameter is on the 
same order of magnitude as the scattering (turbid) medium thickness. 

Application of DSI to the geometry illustrated in Figures 5 and 6 was evaluated numerically 
and experimentally. As shown in figure 6 the scattering medium is a 5 cm-thick slab 
20 containing three identical inclusions located at 5, 25 and 45 mm from the entrance face of the 
slab. The inclusions are at a distance of 4 cm from one another in a direction parallel to the 
face of the slab. The slab is 1 0 cm long (figure 5) with the inclusions spaced 4 cm apart from 
the nearest neighbour inclusion. 
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For the numerical simulations, the turbid medium was assumed to be purely scattering (n, = 
1 mm '^ fia = 0 nun *). The inclusions were asstuned to be point-like, i.e. to have a 
vanishingly small volume. They were also assumed to differ from the scattering medium only 
in their absorption coefficient, such that (volume)i„ciusion multiplied by (^Jinchiaon 5 mm^. 



diffusion model (see M. S. Patterson, B. Chance, and B. C. Wilson, "Time-resolved 
reflectance and transmittance for the noninvasive measurement of tissue optical properties", 
Appl. Opt. 28, 2331-2336 (1989)), adapted to account for the refractive index mismatch 
between the turbid medium and the surrounding environment (see D. Contini, F. Martelli, and 

10 G. Zaccanti, "Photon migration through a turbid slab described by a model based on diffusion 
approximation. I. Theory", Appl. Opt. 36, 4587-4599 (1997)), was used to describe light 
propagation through the simulated medium. The effect of the placement of an inclusion into 
the medium was modelled using a perturbation analysis ( see M. Morin, S. Chatigny, A. 
Mailloux, Y. Painchaud, and P. Beaudry, "Tinie-domain perturbation analysis of a scattering 

15 slab", in Optical Tomography and Spectroscopy of Tissue III, B. Chance, R. R. Alfano, and 
B. J. Tromberg, eds., Proc. SPIE 3597, 67-78 (1 999)). The calculations assumed that 
infinitely short pulses of light were launched into the scattering medium. The transmitted 
light intensity calculation was limited to photons exiting the medium 2 ns after they entered 
the slab. The spatial intensity profiles P^, Pb and Pc thus obtained with each detection 
.20 configuration were calculated for a transverse scan along a line joining the inclusions; P^ is 

profile for configuration 1, Pg is profile for configuration 2 and Pc is profile for configuration 
3. These profiles were then normalized and combined to enhance depth discrimination, as 
discussed above. The results of these combinations are shown in Figure 7 after their 



5 



The surface areas of both the large-area detector and source were set at 400 mm^. The 
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normalization to unity. Figure 7 is a graphic illustration of Spatial profile combinations used 
to enhance longitudinal discrimination. The origin of the abscissa refers to the position of the 
middle inclusion (see Figure S: Front View). The value of Xscm refers to the displacement with 
respect to this origin along the front surface of the medium (see Figure 3: Front View). 

5 refers to the spatial profiles calculated with scanning system configurations shown in figure 2, 
3 and 4 (i.e. namely, configurations 1, 2 and 3). In Figure 7 for the profile of minus Pb is 
illustrated by the graph line which contains the symbols — o— (designated by the reference 
number 30 ); the profile of P^ minus Pc is illustrated by the graph line which contains the 
symbols — (designated by the reference niunber 32 ); the profile of Pb plus Pc minus Pa 

10 is illustrated by the graph line which does not contain any symbol (designated by the 
reference number 34 ). 

Referring to figure 7» as expected, when the profile for configuration 3 (Pc) is subtracted from 
that of configuration 1 (P^), only the signal from the inclusion near the input surface remains. 
On the other hand, subtracting the profile for configuration 2 (Pb) from that of configuration 1 
15 (Pa) nearly eliminates the signal change caused by the inclusion close to the input surface. 
Figure 7 also shows the effect when the profiles for configurations 2 (Pb) and 3 (Pc) are 
added together to get a reisultant sum and the profile for configuration I (Pa) is subtracted 
from the resultant sum. These results clearly show how DSI permits to determine whether an 
inclusion is located close to either surface of the scattering medium. 

20 

The geometry in Figure 5 was investigated experimentally using a phantom consisting of 
three inclusions placed into a 50 mm-thick cell containing a scattering solution of polystyrene 
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microspheres. The inclusions were purely absorbing, 5 mm cylinders of black Delrin and 
positioned according to the geometry illustrated in Figure 5. The scattering solution was non- 
absorbing with ^ = 1 mm *. The laser light was carried close to the scattering medium using 
an input optical fiber. Similarly, the hght transmitted through the mediinn was collected with 
5 an output optical fiber that carried the light to the detector. Using the divergence of light 

emerging from a fiber, a large area source was simulated by simply moving the input fiber ifar 
fix)m the medium. Similarly, a large area detector was simulated by moving the output fiber 
far from the medium. Moving the input fiber away firom the phantom such that photons were 
injected over a surface of 3.1 cm^ functioned as the large area source. The large area detector 
1 0 was likewise achieved by moving the output fiber away fi-om the phantom such that photons 
were collected over a surface of 3.8 cm^. 

The three example optical fiber arrangements are illustrated in Figure 8, 9 and 10. Figure 8 
schematically illustrates ah optical fiber arrangement for a point-soiu-ce/point detector 
configuration, i.e. the optical fibers are configured and spaced so as to provide a scanning 

15 region or surface area and a detection region or surface area which have a small size and are 
more or less the same size. Figure 9 schematically illustrates an optical fiber arrangement 
for a point-source/large-area detector configuration i.e. the optical fibers 40 and 42 are 
configured and spaced so as to provide a scaiming region or surface area and a detection 
region or surface area which are not the same size i.e. the detection area is larger than in 

20 figure 8 - only the output fibers 42 are moved away from the phantom. Figure 10 

schematically illustrates an optical fiber arrangement for a large-area source/point-detector 
configuration i.e. the optical fibers are configured and spaced so as to provide a scanning 
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region or surface area and a detection region or surface area which are not the same size i.e. 
the source area is larger than in figure 8 - only the output fibers 40 are moved away from the 
phantom. Although the figures show that the size of the input and output areas or regions 
may be varied by displacing the fibres 40 and 42 an other type of mechanism may of course 

S be used to effect the change in area. The phantom was scanned using each of these 
configurations using a near-infirared imaging system as described by Y. Painchaud^ A. 
Mailloux, M. Morin» S. Verreault and P. Beaudiy, **Time*domain optical imaging: 
discrimination between scattering and absorption", AppL Opt. 38, 3686-3693 (1999)). The 
scan was done by moving synchronously both the input and output fibers 40 and 42 

10 respectively along the surface of the medium in two dimensions in order to sequentially 

illuminate the entire surface. The detected signal was digitized and sent to a computer that 
was used to construct the images based of the series of measiu-ements performed during a 
scan, each measurement corresponding to one pixel of the image. 

Figure 11, 12 , 13 and 14 show the resultant projection images from scans of a liquid 
1 5 phantom. The horizontal axis corresponds to the transversal placement of the inclusions. 

These images were normalized and combined as described earlier for the simulations. The 
point detector/point source configuration resulted in the image in Figure 14. Figures 1 1 and 
12 show images demonstrating the success of the DSI method in extracting inclusions located 
near the surfaces. A comparison of the images in figures 13 and 14 shows that the inclusion 
20 placed in the central plane of the phantom was highlighted relative to those near the surface 
using DSI. 
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As may be appreciated from the above, the Dual Spatial Integration technique (DSI) of the 
present invention is proposed as a tool to aid in the visualization of structures located in the 
central plane of a turbid medium. DSI represents an improvement with respect to the prior art 
in that one gains depth information content usually obtained with multi-point detection 
5 geometry without its concomitant time costs. This is achieved by using different source- 
detector configurations to section the medium, namely as follows: 

1 . The scattering medium is scanned using a point source and point detector (i.e* 
configuration 1). 

2. An image (I a) is generated based on the detected light using configuration 1 . As an 
10 example, the detection of light can be done in a time*resolved manner in which only the first 

arrival photons are considered. The image can also be generated using indirect characteristics . 
of the detected light as in the FIDM method described elsewhere^ 

3. The scattering medium is scanned using a point source and a large area detector 
(configuration 2). 

15 4. A second image (Ib) is generated using the same method as in 2 but with 

configuration 2. 

5. The scattering medium is scanned using a large area source and a point detector 
(configuration 3). 

6. A third image (I^) is generated using the same method as in 2 but using the 
20 configurations. 
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7. Image Ib is subtracted fix>m image 1^ to enhance structures located close to the 
output surface of the scattering medium. 

8. Image Ic is subtracted from image 1^ to enhance structures located close to the input 
surface of the scattering medium. 

5 9. Image Ia is subtracted from the sum of images Ib and Ic to produce an image 

wherein the centrally-located stmctures are enhanced relative to those close to the siurface. 

Note that other combinations of images 1^, Ib wd le or other data processing involving 1^, Ib 
and Ic may also result in obtaining relevant information. 

In the summary given above, a point source refers to the injection of light in a medium over 
10 a small area; a large area source refers to the injection of light in the medium over a large 
area; a point detector refers to the detection of light that emerges from the medium over a 
small area; a large area detector refers to the detection of light that emerges from the medium 
over a large area. Note again that "small area" refers to a siurface whose diameter is small 
compared to the medium thickness while "large area" refers to a surface whose diameter is on 
15 the same order of magnitude as the medium thickness. The appropriate input and output area 
sizes may be determined empirically for any given turbid mediiun to be subjected to analysis. 
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We claim: 

1 . In a method for scanning a turbid medium for generating an image thereof for the 
detection of one or more anomalies contained within the turbid medium, said turbid 
medium having a first face and an opposite second face, wherein said turbid medium 
5 is scanned by displacing an optical signal source and a corresponding optical detector 

from one respective spatial location to another about the turbid medium, each spatial 
location being associated with a corresponding input region on said first face and a 
corresponding output region on the opposite second face, said optical signal source 
directing an optical signal to each of said input regions, said optical detector in 
10 response to optical signals detected from each of said output regions sending 

corresponding detector output signal data to an image processing means for generating 
a set of image data of the turbid medium, the improvement comprising 

generating a primary set of image data derived by scanning said turbid medium using 
for each input region a first input region of predetermined size and using for each output 
15 region a first output region of predetermined size, and 

generating a secondary set of image data derived by scanning said turbid medium 
using for each input region said first input region of predetermined size and using for each 
output region a second output region of predetermined size greater than that of said first 
output region. 

20 2. A method as defined in claim 1 comprising generating a composite set of image data 
by subtracting the secondary set of image data from the primary set of image data. 
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3. In a method for scanning a turbid medium for generating an image thereof for the 
detection of one or more anomalies contained within the turbid medium, said turbid 
medium having a first face and an opposite second face, wherein said tuibid medium 
is scanned by displacing an optical signal source and a corresponding optical detector 
S from one respective spatial location to another about the turbid medium, each spatial 

location being associated with a corresponding input region on said first face and a 
corresponding output region on the opposite second face, said optical signal source 
directing an optical signal to each of said input regions, said optical detector in 
response to optical output signals detected from each of said output regions sending 
10 corresponding detector output signal data to an image processing means for generating 

a set of image data of the turbid medium, the improvement comprising 

generating a primary set of image data derived by scaiming said turbid medium using 
for each input region a first input region of predetermined size and using for each output 
region a first output region of predetermined size, and 

IS generating a further set of image data derived by scanning said turbid medium using 

for each input region a second input region of predetermined size greater than that of said first 
input region and using for each output region said first output region of predetermined size. 

•4 A method as defined in claim 3 comprising generating a composite set of image data 
by subtracting the further set of image data firom the primary set of image data. 

20 5. In a method for scanning a turbid medium for generating an image thereof for the 
detection of one or more anomalies contained within the turbid medium, said turbid 
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medium having a first face and an opposite second face, wherein said turbid medium 
is scanned by displacing an optical signal source and a corresponding optical detector 
from one respective spatial location to another about the turbid medium, each spatial 
location being associated with a corresponding input region on said first face and a 
5 corresponding output region on the opposite second face, said optical signal source 

directing an optical signal to each of said input regions, said optical detector in 
response to optical output signals detected firom each of said output regions sending 
corresponding detector output signal data to an image processing means for generating 
a set of image data of the turbid medium, the improvement comprising 

10 generating a primary set of image data I a derived by scanning said turbid medium 

using for each input region a first input region of predetermined size and using for each 
output region a first output region of predetermined size, 

generating a secondary set of image data Ib derived by scanning said turbid medium 
using for each input region said first input region of predetermined size and using for each 
1 5 output region a second output region of predetermined size greater than that of said first 
output region and 

generating a tertiary set of image data 1^ derived by scanning said turbid medium 
using for each input region a second input region of predetermined size greater than that of 
said first input region and using for each output region said first output region of 
20 predetermined size. 

6. A method as defined in claim 5 comprising generating a composite set of image data 
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5 



8. 



10 9. 



in accordance with the equation 1b + Ic - 

A method as defined in claim 1 comprising combining said primaxy set of image data 
and said secondary set of image data so as to generate a product set of image data 
wherein the similarities and/or differences between the said primary set of image data 
and said secondary set of image data are highlighted. 

A method as defined in claim 3 comprising combining said primary set of image data 
and said further set of image data so as to generate a product set of image data 
wherein the similarities and/or differences between the said primary set of image data 
and said further set of image data are highlighted. 

A method as defined in claim S comprising combining said primary set of image data 
said secondary set of image data and said tertiary set of image data so as to generate a 
product set of image data wherein the similarities and/or differences between said 
primary set of image data, said secondary set of image data and said tertiary set of 
image data are highlighted. 
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